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^ We discuss possible origins for discrepancies observed in the radiated energies in head-on 

, collisions of non-spinning binaries starting from Brill-Lindquist and superposed Kerr- 

Schild data. For this purpose, we discuss the impact of different choices of gauge param- 
eters and a small initial boost of the black holes. 

1. Introduction 



The area of numerical relativity has made dramatic progress in the last two years'^'^'^ 
and numerical simulations of black hole binaries performed by various groups have 
resulted in a wealth of literature on astrophysical topics and those related to gravi- 
tational wave data analysis. At the same time, laser-interferometric GW detectors, 
^1 LIGO, GEO600, TAMA, VIRGO, have started collecting data at design sensitivity. 

The area of GW physics has thus entered a very exciting stage with vast potential 
for astrophysics, our understanding of the early universe and fundamental physics. 
^ ' From the viewpoint of numerical relativity, though, a number of important ques- 

, tions still remain to be addressed. These largely concern the accuracy of the pro- 

I duced waveforms, the dependency of the results on numerical techniques, their de- 

' tailed matching with results predicted by approximation theories as well as the 

mass production of waveforms covering the complete parameter space for use in 
' GW observations and parameter estimation. 

. The purpose of this study is to address the dependency of the numerical results 

on the choice of black-hole binary initial data. In contrast to spacetimes containing 
single stationary black holes, there exist no uniqueness theorems guaranteeing that 
two data sets for binary black holes using different data types represent the same 



o 
o 



X 



5h ' physical configuration. Indeed, such data sets are known to generally differ in the 

^ amount of gravitational radiation inherent in the initial data. 

The dependency on initial data parameters (though not data type) has been 
studied in the case of binary black hole coalescence in Refs. 2,4,6. Using the moving 
puncture technique, the merger waveforms are found to agree well for different 
initial separations and algorithms to produce quasi-circular initial configurations. A 
comparison of GWs produced in the evolution of Cook-Pfeiffer and puncture data 
using different evolution techniques has been presented in Ref. 1 and shown good 
agreement. There remains a difficulty in the identification of free initial parameters 
in this case, however (cf. the non-vanishing spin in the Cook-Pfeiffer data set in 
this comparison). This identification of parameters represents a simpler and cleaner 
task in the case of head-on collisions of non-spinning black holes which has been 
studied in Ref. 11. That study observed systematically larger amplitudes by about 
10 % in the merger waveform resulting from Kerr-Schild data compared with those 
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of Brill-Lindquist data. Here we investigate two possible causes for this discrepancy: 
the dependency of the results on the gauge trajectories in the case of Kerr-Schild 
data and the impact of deviations from time symmetry of the initial data. 

2. Results 

The simulations presented in this work have been obtained with the Lean code''''^"'^ 
which uses the BSSN formulation of the Einstein equations together with the moving 
puncture approach. '^'^ It is based on the Cactus^ computational toolkit and the 
Carpet-'^° mesh-refinement package. For a detailed description of the code as well 
as the construction of initial data we refer the reader to Ref. 11. 

We first discuss the gauge trajectories used in Ref. 11 for the Kerr-Schild data. 
There, algebraic gauge conditions are constructed which require trajectories for the 
(approximate) black hole positions (see 11 for details). These are prescribed as 
polynomials ±x''{t) = Xq + v^t + + + g't"'/24 which are smoothly (up 

to fourth derivatives) matched to the static function a;*(t) = in a time interval 
ti < t < t2- Here Xq, Dq, Cq, ^q, ti and t2 are free parameters which need to be 
chosen carefully to avoid numerical instabilities. In Table [T] we list the values for 
each (Kerr-Schild) model. In order to assess the impact of the particular choice of 
these parameters, we have evolved the initial data of model 1 with alternative gauge 
parameters as listed in the second row of the table. This alternative gauge trajectory 
is motivated by the initial coordinate velocity v — —0.08 of the central position of 
the apparent horizon as measured using Thornburg's AHFinderDirect.^'^'^'^ 

The resulting waveforms are shown in the left panel of Fig.[TJ Both the waveforms 
and the radiated energies thus obtained for model 1 show excellent agreement. The 
differences in radiated energy are about 1.5 % and thus substantially smaller than 
the discrepancies between ii^KS and -Ebl- 

Second, we assess the impact of deviations from exact time symmetry of the 
initial superposed Kerr-Schild data. These deviations manifest themselves in a small 
but non- vanishing initial coordinate velocity of the superposed Kerr-Schild holes as 
measured by the central position of the apparent horizon. For the case of model 2 

Table 1. Parameters for the black hole models. The gauge parameters in columns 5 to 9 are 
only used for the Kerr-Schild simulations. There we also use = —0.000278, —0.000165 and 
—0.000104 respectively for models 2a, 3 and 4. Here M is the ADM mass of the system. 
The two rightmost columns list the energy radiated in the ^ = 2, m = mode for both data 
types, ignoring contributions due to the spurious initial burst. 
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0.067 
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we have measured this velocity to be w = 0.067. In order to estimate what impact 
such an initial velocity has on the resulting waveforms, we have applied an initial 
linear momentum pz = m v to the Brill-Lindquist version of this model, where m is 
the irreducible mass of a single hole. The resulting waveform is compared with its 
non-boosted counterpart in the right panel of Fig. [TJ Again, the wave amplitudes 
show good agreement, as do the resulting values for the radiated energy in Table [TJ 
In summary, we find the observed differences in radiated energy resulting from 
modifications of the gauge trajectories and a possible initial boost of the black 
holes to be of the order of 1 % and thus substantially below the differences of about 
> 20 % observed in Ref. 11 between the two types of initial data. 
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